function vf = forwardSpeed(x, vmax)
  t = 0.27;
  k = 0.19;
  vf = k * x;
  vf = vf .* (vf > k * t);
  vf = vf .* (vf < vmax) + vmax .* (vf >= vmax);
end
